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ABSTRACT: We consider twist-3 operators in the s((2) sector of M = 4 SYM built out 
of three scalar fields with derivatives. We extract from the Bethe Ansatz equations of this 
sector the exact lowest anomalous dimension 7(5) of scaling fields for several values of the 
operator spin s. We propose compact closed expressions for the spin dependence of j(s) 
up to the four loop level and show that they obey a simple new twist-3 transcendentality 
principle. As a check, we reproduce the four loop universal cusp anomalous dimension 
governing the logarithmic large spin limit of 7(5). 
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1. Introduction 

Integrability in QCD is an old intriguing issue (see [|l]] for a recent review). Certainly, QCD 
is not an integrable quantum field theory. However, in the planar multicolor 't Hooft limit, 
new special features emerge due to the simplified dynamics. For instance, it is known 
that the one-loop renormalization group flow of certain sets of composite operators is 
associated with integrable XXX lattice Hamiltonians with sl(2) spin symmetry The 
evolution time in the integrable model is the (logarithm of the) renormalization scale at 
which we define the composite operators. The integrable Hamiltonian is identified with 
the quantum dilatation operator. Its eigenvalues are the anomalous dimensions of the 
scaling fields. Several explicit examples can be found in [[| ||, |j| |6[ |7| |8[ |10| ]. 

Inspired by the parallel developments of AdS/ CFT duality j^], the integrability prop- 
erties of large N gauge theories have been deeply investigated in supersymmetric Yang- 
Mills theories. In this paper, we shall work in the context of the maximally supersymmet- 
ric M = 4 super Yang-Mills theory which is UV finite and superconformally invariant at 



-1- 



the quantum level. In the analysis, we shall always understand the planar limit. Due to 
the possible interdisciplinary interest of our investigation, we briefly overview the main 
logical developments of integrability in M = 4 SYM. 

The first positive results are described in the seminal paper [|Tl|]. The one loop dilata- 
tion operator is computed in a closed so (6) sector of scalar operators. It can be identified 
with the Hamiltonian of an integrable lattice model. Soon, the integrable structure could 
be extended to the full set of psu(2, 2|4) operators. Also, it appeared clear that integrability 
could work beyond one loop, with the coupling constant being a deformation parameter 
for the superconformal algebra representation [12, [jj, 14, 15]. 

These investigations suggested that M = 4 SYM could admit an all loop solution for 
the renormalization flow of its composite operators. The reason for this peculiar internal 
integrability is the strong constraint imposed by the unbroken superconformal symme- 
try Also, in agreement with AdS/ CFT duality, the symmetry is further enhanced at tree 
level to the higher spin algebra fjs(2, 2|4). Its multiplets obey intricate recombination rules 
under the interacting subalgebra psu(2, 2|4) [|T6|]. 

The analysis of the dilatation operators was further extended in [|T7|], computing the 
three loop dilatation operator in the 5u(2|3) subsector of the theory. This remarkable 
achievement has been obtained algebraically by exploiting the superconformal invariance 
instead of direct Feynman diagram calculations. The su(2 1 3) sector includes closed smaller 
subsector s like the bosonic su(2) and the fermionic su(l|l), but not the non-compact s((2) 
containing the twist operators. 

As is common in the context of integrable models, the basic object is not the Hamil- 
tonian, but the 5-matrix for elementary excitations. Integrability means that scattering 
processes are elastic and factorized. To a large extent, the two body scattering matrix al- 
lows to reconstruct the full dynamics. Once we know the 5-matrix, we can write down 
Bethe Ansatz equations that determine the allowed quantum numbers and momenta of 
the elementary excitations. From the solutions of these equations we extract the spectrum 
of the integrable Hamiltonian, i.e. the anomalous dimensions in the gauge theory. 

Therefore, as soon as integrability was recognized as an essential feature of M = 4 
SYM, a new investigation route started with the aim of computing higher order Bethe 
Ansatz equations bypassing the explicit construction of the dilatation operator. For in- 
stance, a remarkable result was obtained in [ |T8| ] where the three loop Bethe Ansatz equa- 
tions were determined in the su(2) C su(2|3) bosonic sector. 

Here and in the following, the Bethe Ansatz equations that we shall discuss will al- 
ways be asymptotic. This means that their perturbative expansion predicts the correct 
anomalous dimension of the operator O up to a finite loop order 0(g 2L ). The number L 
depends in a controlled way on the number of elementary fields in O. The dilatation op- 
erator, seen as a spin chain Hamiltonian, contains long-range interactions up to a distance 
growing with the perturbative order. The number of fields in O determines the length of 
the spin chain and, beyond a certain order, wrapping problems appear invalidating the 
Bethe Ansatz equations. 
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Due the simplicity of the su(2) sector, it was possible to compute the three loop dilata- 
tion operator and its Bethe Ansatz equations at the five loops level []19|]. This remarkable 
result was also confirmed for many operators by standard Feynman diagram calculations 
at three loops p0|]. The Bethe Ansatz and S-matrix approach culminated in the work 
A clever combination of gauge /string theory arguments allowed to conjecture the three 
loop Bethe Ansatz equations for the other rank-1 closed subsectors su(l| 1) and sl(2) with- 
out computing the dilatation operator. 



From our point of view, the paper ply is crucial because it started the investigation of 
the non-compact sl(2) sector. As a check of the proposed Bethe Ansatz equations, a few 
explicit two loops calculation have been performed in [Q] using superspace perturbation 
theory. Also, the two loop dilatation operator has been constructed algebraically in the 
5u(l, 1|2) D s[(2) sector p3|]. In all cases, the Bethe equations of [ pip were fully confirmed. 

The three loop 5-matrix in the s[(2) sector is an important bridge toward QCD. Based 
on three loop QCD calculations [[24p and inspired by one and two loop results in JV = 4 
SYM p5| , ^6p, Kotikov, Lipatov, Onishchenko and Velizhanin (KLOV) conjectured a three 
loop prediction for the anomalous dimension of JV = 4 twist-2 superconformal opera- 
tors at generic spin in [g7|]. The prediction is based on what is now called the maximum 
transcendentality principle. 

This prediction is perfectly matched by the perturbative expansion of the sl(2) Bethe 
Ansatz equations. Notice that at twist-2 all conformal operators fall in a single supermul- 
tiplet (see [M]]). With more complicated operators this is not necessarily true. For 



instance, the paper piQ studies the two loop dilatation in M = 1,2,4 SYM for Wilson op- 
erators with 3 quark /gaugino fields and derivatives. This is not the s((2) sector that we 
are discussing and that is built with the holomorphic scalar fields of N = 4 SYM. 

As a technical remark, the agreement between [ pT| ] and the KLOV prediction is some- 
what beyond the regime of applicability of the Bethe equations as it would follow from 
wrapping considerations. This subtle point is related to superconformal symmetry as ex- 



plained in the later works [ |32| , |33| ]. Wrapping problems are expected to appear at L + 2 



loop order for twist- L operators. 

Soon, it was realized that the wrapping barrier could be overcome by considering 
the large spin limit of anomalous dimensions. In this limit, the leading contribution to 
the smallest anomalous dimension scales logarithmically with the spin s and is universal 



with respect to twist |3£|. The coefficient of log s is a non-trivial function of the coupling, 



the cusp anomalous dimension (a.k.a. scaling function) [p^pq, 



Since wrapping problems do not affect the calculation of the scaling function, a pos- 
sibility opened for its all-loop calculation from the Bethe Ansatz equations. The first at- 



tempt in this direction is described in u38fl . In that work, the sl(2) Bethe equations are put 
on a more solid basis, by a field theoretical calculation of the two loop dilatation operator. 
This calculation confirmed that it is diagonalized by the Bethe Ansatz equations of [ffifl]. 
Then, an integral equation was proposed to determine at all loops the scaling function, 
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the Eden-Staudacher (ES) equation. As a matter of fact, the four loop ES prediction vio- 
lates the maximum transcendentality principle. This anomalous fact was soon related to 
a missing crucial piece in the Bethe Ansatz analysis, namely the dressing factor [|39|]. It is 
an Abelian phase which is not constrained by superconformal symmetry and enters the 
perturbative expansion of anomalous dimensions precisely at four loops. 

Although, at first, the appearance of the dressing factor is rather discouraging, AdS/ CFT 
duality comes to rescue. The M = 4 SYM theory has a AdSs x S 5 string counterpart which 
is also classically integrable (and to some extent also at the quantum level) and where an 
all-order strong coupling expansion for the dressing is known. This is made possible by 
a combination of semiclassical string calculations and string integrability considerations, 
like crossing symmetry [||, |T|, ||, |3[ ||, ||]. 



By an impressive insight, Beisert, Eden and Staudacher proposed in pq ] an all-order 
weak-coupling expansion of the dressing phase. The modified integral equation for the 
scaling function now predicted a four loop contribution in analytical form in agreement 
with the KLOV maximum transcendentality principle. In a synchronized but independent 
fashion, a remarkable four-loop Feynman diagram calculation appeared p7[ ] providing a 
numerical expression for this contribution in full agreement with the result of p6|]. 

More recently, there have been additional important developments aimed at a deeper 
understanding of the dressing factor Ji^ , |4l| 50, 51]. In this paper, we shall dwell on the 
results of [f§6|] to compute certain 4 loop contributions. 

A complementary approach to integrable systems (closely related to the Bethe Ansatz) 
is based on the Baxter Q-operator [p^]. In this framework, the two loop dilatation oper- 
ator in M = 2, 4 SYM for Wilson operators with scalars and derivatives (our st(2) sector) 
as well as the two loop Baxter operator are computed in []53p. The all-loop asymptotic 
generalization of the Baxter equation in the same sector also appeared in [B3]. All-loop 
extensions to the larger sl(2| 1) sector have been recently published in [ |55| , |55| ]. 

In this paper, we keep working in the s((2) sector of TV = 4 SYM and study twist-3 
operators extracting from the Bethe Ansatz equations the (lowest) anomalous dimension. 
For these operators, the Bethe Ansatz equations plus dressing are expected to be reliable 
up to 4 loops, even at finite spin. We shall show that, for even spin, the minimal anoma- 
lous dimension is associated with an unpaired operator and the non-dressing part of it 
has rational contributions up to 4 loops. 

We shall be able to provide closed expressions for the various perturbative contribu- 
tions at finite spin s, including the dressing part at four loops. As a non trivial check, we 
shall recover the universal four loop cusp anomalous dimension. 

The detailed plan of the paper is as follows. Sec. (|l]) is devoted to a brief overview of 
the main integrability facts in M = 4 SYM with some discussion of the links with QCD 
integrability. Sec. @ describes the relevant Bethe equations. Sec. @) is devoted to the 
detailed presentation of our analysis on the twist-3 operators. A self-contained Appendix 
contains some technical information about nested harmonic sums, which are the basic 
element entering our proposed closed formulae for anomalous dimensions at finite spin. 
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2. The st(2) sector of M = 4 SYM 



2.1 The Bethe Ansatz equations 

The sl(2) sector of planar J\f = 4 SYM contains single trace states which are linear combi- 
nations of the basic operators 

Tr {(V s ' Z) ■■■ (V s l Z)}, Sl + ... + s L = s, (2.1) 

where Z is one of the three complex scalar fields and V is a light-cone covariant derivative. 
The numbers {s{\ are non-negative integers and s is the total spin. The number L of Z 
fields is the twist of the operator, i.e. the classical dimension minus the spin. The subsector 
of states with fixed spin and twist is perturbatively closed under renormalization mixing. 

At one-loop, the dilatation operator in this sector maps to an integrable spin chain 
with L spins transforming according to the s = — 1 infinite dimensional sl(2) represen- 
tation Jl3|]. Beyond one loop, the work [[n]] proposed asymptotic all-order Bethe Ansatz 
equations. As a check, a few explicit two loops calculation have been performed in [22] 



using superspace perturbation theory. Also, the two loop dilatation operator has been 
constructed algebraically in the su(l, 1|2) D s[(2) sector |23|]. It has also been evaluated in 



the s 1(2) by Feynman diagram calculations in [p8[]. In all cases, the Bethe equations of [ |2l| ] 
were confirmed. Wrapping problems are delayed by supersymmetry and appear at L + 2 
loop order for twist-L operators [ 02] , [33| ] (also, N. Beisert and M. Staudacher, private com- 
munication). 

The anomalous dimensions of scaling combinations of states Eq. ( |2.1[ ) are the eigen- 
values 7l(s; g) of the dilatation operator /integrable Hamiltonian where our definition for 
the planar coupling is 

9 2 = ( 2 - 2 ) 

where N is the number of colors. The coupling g is kept fixed as N — > 00. These anoma- 
lous dimensions 7l(s; g) are obtained by solving perturbatively the Bethe Ansatz equa- 
tions, provided we are in the wrapping-free cases. 

The Bethe Ansatz equations determine s real Bethe roots {v,k}i<k<s and read 

, 9 2 



=n i: x j- 2x f , 4=*(u k ±-), (2.3) 



, X k J j^k X k ~ X j 1 9 " 2 



2 x k Xj 



where we have defined the maps 



2g 



g 2 



x{u) = -\ l + \\l-^ , U {x) = x + ± z . (2.4) 



2x 



The relevant solutions of Eq. ( JZ.3| ) are those obeying the constraint Y\k=ii x k / x k ) = ^ 
needed to project onto cyclic states associated with single-trace operators in the gauge 
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theory. Given a solution of Bethe Ansatz equations, we obtain the anomalous dimension 
from the formula 

7,(^)= 5 2 tf^-^)- (2-5) 

i=i \ i % I 

Often, the Bethe Ansatz equations are presented in a different, equivalent form, more 
suitable for our later discussion. We introduce the Bethe momenta {pi}i<i< s which are 
related to the Bethe roots by the transformations 

p(u) = -ilog^P-, (2.6) 
X (u) 



U(P) = ^cot^/l + 8 5 W^ 



(2.7) 

The x^ combinations appearing in the Bethe equations have the following expression in 
terms of the p;, 

_j_ ip . 

I± < I " = ^( 1+ V 1+8! ' 2sin2 f)- (Z8) 

Finally, the anomalous dimension reads 



ids- g) = ( \/l + 8 5 2 sin 2 ^ - 1) . (2.9) 



Notice that in this form, the anomalous dimension is obtained as a sum of single quasi- 
particle energies with a definite dispersion relation. 

The Bethe equations admit several solutions that describe different scaling operators. 
The different solutions are selected by taking the logarithm of the Bethe equations and 
choosing a determination for the branch. This is known as choosing the mode /Bethe 
numbers of the state. In the following analysis we shall study the minimal anomalous 
dimension, i.e. the ground state of the integrable spin chain, for reasons that shall be- 
come clear in the discussion. The solution describing the ground state has known mode 
numbers at any twist [f38|]. 

Starting from the 4-loop level, the Bethe equations Eqs. (|Z3|) must be modified includ- 
ing a universal Abelian dressing phase. In other words, beyond three loops, the correct 
equations are 

9 2 

. . . r. . 1 - 



[£*.) = rr^^ — x -^e 2 ^, (2.io) 
\x k J ^ k H-x 31 _JL_ 

x k x j 

with tfkj = 0(g 6 ). The general perturbative expansion of the dressing phase takes the 
form 

*H = £ E E T [9r(w)9r + l+*(Pi) - (* « J)} > (2-H) 

r>2i/>0/j>i/ V / 
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where the higher order charges q r (p) have the expression [ [T^ ] 



2sin(^V) /yi + 8 5 2 sin 2 f - 1 
r — 1 \ 2 g 2 sin ■ 



r-1 



(31 

with the non trivial coefficient being 3 7^ 0. 



The proposed formula for the coefficients p£ h is given in [ }46| ] (see also |57l £58|] for 



related explorations) and reads 



Jr+u+n) _ , lY+ » + i (r-l)(r + 2u) ( 2^ + 1 \ /2 M + l\ 



+ 1), (2-13) 



and zero if /i — r — z/ + 1 < 0. In particular, for the leading order 4-loop correction we have 

2# kj = (qi°\p k ) - q { °\ Pj ) q { °\p k )) g 6 + 0{g% (2.14) 

g( 0) (p) = 4sin 2 |, (2.15) 

(p) = 4 sin 2 I sinp, (2.16) 



and 



/3g =4/3, /3 = Ca. (2-17) 



2.2 Extracting the perturbative anomalous dimensions, the twist-2 case 

Let us fix the twist. For a given spin, we can start by solving numerically the one-loop 
Bethe Ansatz equations (see later for more analytical information on this step). Then, it is 
straightforward to work out the perturbative expansion of the all-order equations. In this 
way we obtain numerical values for the various loop contributions to 



7L(s;5) = E7i n W"- (2-18) 
n>l 



If these coefficients are rational, we can identify them unambiguously by working with a 
very large number of digits. 

Actually, in many integrable systems, there is an alternative approach based on the Q 
Baxter operator [|52|]. One build a polynomial Q(u) whose roots are the Bethe roots. This 
polynomial satisfies a recurrence equation that determines it completely, once the quan- 
tum numbers of the desired state are chosen. From the (analytical) knowledge of Q(u), 
it is possible to obtain the Hamiltonian eigenvalues without resorting to finite precision 
numerical calculations. 

At twist-2, the one loop Baxter operator is known in closed form and the Baxter equa- 
tion is also known up to three loops. The Bethe Ansatz and Baxter method of course agree 
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as we checked up to spin s = 68. The resulting anomalous dimensions agrees with the 
KLOV three loop prediction 

7?] =4 Si, (2.19) 



72 



(2) 



-4 (s 3 + S_ 3 - 2 S_2,i + 2 Si (S 2 + S_ 2 ) 



7$ = -8(2 5_ 3 S 2 - S 5 - 2 5_ 2 S 3 - 3 S_ 5 + 24 5_ 2 ,i,i,i 

+ 6 (S_4_i + S_3_2 + S_2,3) — 12 (S-3^1 + S_2,l,2 + S_2,2,l) 

- (S 2 + 2 S 2 ) (3 S_ 3 + S 3 - 2 S_ 2 ,i) - Si (8 S_ 4 + S 2 2 
+ 4 5 2 S_ 2 + 2 S 2 + 3 S 4 - 12 5_ 3 ,i - 10 S_ 2 , 2 + 16 S_ 2 ,i,i 



In the above expressions, all the S functions are harmonic functions (see Appendix A) 
evaluated at the argument s, the spin 

S a = S a (s). (2.20) 

This check is not surprising and extend the checks already performed in [p]]. 
3. The sl(2) sector at twist-3 

The same approach can be applied to the twist-3 case. Here, the one loop anomalous 
dimension 73 (s) is known from the analysis of |14] (see also some conjectures in flT6[l). 



In particular, for even spin, the ground state of the integrable Hamiltonian is an unpaired 
state. 

At twist-3, the Bethe equations are believed to be reliable up to four loops. This is 
beyond the range of applicability of the known Baxter equation that stops at three loops. 
So, we have computed 73 1,2 ' 3 ^(s) by both methods up to s = 68 and 73^(5) in the same 
range by the Bethe equations only. 

(n) 

This means that we have computed 73 '(s) for n = 1, 2, 3, 4 and 5 < 70 in analytical 
form as rational numbers plus, at four loop, a transcendental dressing contribution. The 
limit s < 70 is fixed by computer limitations that appear when the result is sought in 
analytical form. At s = 68 and four loops we need several thousands of digits. Of course, 
if a numerical result is enough, it is possible to push quite further the calculation as we 
shall discuss later. 

In the following sections, we shall first present new exact results for the twist-3 Baxter 
operator at one loop, as well as all the details of the Baxter three loop calculation. 

3.1 The one-loop Baxter operator at twist-3 



We follow the notation of [p4[]. The one-loop Baxter equation at twist-3 reads 



u + -\ Q 3 (u + i)+ lu - -) Q 3 {u-z) = t 3 {u)Q 3 {u), (3.1) 
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where 

h(u) = 2u 3 + q 2 u + q 3 , (3.2) 

•=-H)NK (3 - 3) 

Here, s is the (even) spin. The Baxter operator Q(u) is a polynomial of degree s 

s 

Q 3 (u) = J2 a nU n , (3.4) 

71 = 

whose roots are the Bethe roots. Replacing Q 3 in the Baxter equation we obtain a homo- 
geneous linear problem in the coefficients {a n }. The parameter q 3 is a quantum number. 
It appears in the linear problem as an eigenvalue. For the lowest state, 53 = 0. Also, the 
Baxter polynomial turns out to be even under u — > —u. 

The eigenvector associated with q$ = determines {a n } and hence Q3 up to an irrel- 
evant scaling factor. Given the roots {u n } of the Baxter polynomial 



Q 3 (u) = ■ fl (u - u n ), (3.5) 

n=l 

we can compute the one-loop anomalous dimension by the formula 

* )l,) =i*hr*='«tt)- (3 - 6) 

Notice that the last expression does not require the knowledge of {u n }, but just Qs(u). 
For example, at s = 2 we find 

Q 3 (tt)~(2tt + l)(2tt-l), « = ±^> 7i 1} (2) = 4. (3.7) 
For s = 4 we find 



Q 3 (u)~H-72u 2 +48u 4 , u = ±^|±_l 7^(4) = 6. (3.8) 

Now, with some trial and error and guided by a comment about Wilson polynomials 
in [34], we found the following closed solution to the Baxter equation 

/_£ £ I i _ ■ \ 

Qs^A/^s 2'2 + 1 '2 + ' U, 2 . (3.9) 

V 1,1.1 / 

Then, one obtains 

if'M = (-|) = « (f + l) . ^ ( 1 - 5" 2 + 5 ; 1 J = ■ ■ ■ = 45, (i) . (3.10) 
This result is new and is in perfect agreement with the empirical conjecture of [14]. 
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3.2 The three loop Baxter equation at twist-3 

The Baxter equation at twist-3 valid up to three loops is described in Let Q(u) be the 
Baxter polynomial. For simplicity, we omit the twist label. Q can be loop expanded in the 
form 

Q(u) = QW (u) + g 2 Q (1) (u) + 5 4 Q {2) («) + ■■■■ (3.11) 

The polynomial Q(°\u) has degree s and is = Q3 defined in Eq. (f3l^). The polyno- 
mials Q^ k \u), k = 1, 2 have degree s — 2 and, for the ground state, are also even under 
u — » — u. Hence, they have the general form 



n=0 



The three loop Baxter equation reads 



Q(u + i) + A. 



x ( it + - 



x \ u 

1 2 



Q(u-i) = t(u)Q(u), 



where we have defined again 



u 



x(u) = - ( 1 + 



2 5 2 



It* 



(3.12) 



(3.13) 



(3.14) 



and 



A±(i) = a; 3 exp (logQ(z)V 



(3.15) 



t(u) 



4 a: 2 



(log Q(z))^ =±1 +x (log Q(z)) 



in 

z=±i 



A + (x(u)) A_(x(u)) 



2+E 



Qn(g) 



n=2 (*(«))". 

The charges c? n have a weak-coupling expansion that for the ground state is 



(3.16) 



92(5) = - ( 3 + 

33(5) = 0, 

00 

n=l 

q n (g) = 0, n > 4. 



s + 



(3.17) 
(3.18) 
(3.19) 
(3.20) 



Given the (unique polynomial) solution Q(u) of the Baxter equation, the 3-loop anoma- 
lous dimension is given by the formula 



l3 (s; g) = 2 iL 2 (log Q(u))' u ^ + £ (log Q(u))^ + ^ (log Q(u))^ f + 0( 5 8 ) 



(3.21) 
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where we exploited the parity invariance of Q. Notice that the above formula must be 
further expanded in g 2 since the polynomial Q contains the coupling g. 

We did not find an analytical formula for the higher Baxter polynomials Q^ k \ k = 2,3. 
However, it is straightforward to compute them for each value of the spin s. 
For example, at s = 2 we find 

Q = \^u 2 -l)-\g 2 + Q-g\ (3.22) 
73 (2) = 4 5 2 -6 5 4 + 17 5 6 + --- . (3.23) 



At s = 4, 



1 a 2 
Q = — (48 w 4 - 72 u 2 + 11) - — (108 w 2 - 47) + (3.24) 
48 24 

+f-(36u 2 + 53), 
lb 

73(4) = 6 5 2 -f 5 4 + ^ 5 6 + .... (3.25) 



At s = 6, 



Q = — (320 w 6 - 1520 u 4 + 1292 u 2 - 153) + (3.26) 
320 

a 2 

- — (8400 w 4 - 17928 u 2 +4357) + 
800 v ' 

3o 4 

+ — — (178000 u 4 + 325448 u 2 - 359587), 
64000 v ; ' 

22 - 443 4 303115 fi tn nrn 

^ = Y 9 + ^rm g + — (3 - 27) 

3.3 The proposed four loop anomalous dimension at twist-3 

Using the methods described in the previous sections, we obtained the coefficients 73 ' (s) 
as rational numbers for even s < 70, plus a transcendental term at four loop coming from 
the dressing to be discussed below, separately. 

Is it possible to predict a closed formula for the coefficients as functions of s ? Cer- 
tainly, one needs some inspiration to convert the problem into a well-posed one. We 
propose the following claim guided by several numerical explorations: 

Twist-3 transcendentality principle: The expression 0/73 '(s) is obtained as a sum of 
transcendentality 2n—l (products of) nested harmonic sums with positive indices and argument 
s/2. For n = A the same holds with the dressing contribution being £3 times a transcendentality 4 
combination. 

Counting the number of possible harmonic sums (see Appendix (A)), one sees that 
this principle combined with the s < 70 results fixes the three loop anomalous dimension. 
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One finds the remarkably simple expressions 
7^=4*, 

4 2) = -2(5 3 + 2 5!5 2 ) (3.28) 
7 J 3) = 5 S 5 + 6 S 2 S 3 - 8 5 3 ,i,i + 4 5 4 ,i - 4 S 2 , 3 + 
+5i (4 5| + 2 5 4 + 8 5 3 ,i), 

with all 5 a functions are evaluated at the argument s/2, i.e. half the (even) spin 

S a = S a (~) . (3.29) 



,2. 

The argument 5/2 is not surprising as we know from the one-loop analysis. 
At four loops, the anomalous dimension has two contributions 

7 ^ = 7 | 4 ' no dressing ) + 7 ( 4 - dressin g) . (3.30) 

The contribution from the dressing is included in the Bethe equations according to Eq. (2.12) 
of p6|]. The expression for »y( 4 j dressing ) [ s easily found by assuming that /3 has transcenden- 
tality 3. The current Ansatz is j3 = £ 3 . Applying our principle, one fixes 

^dressing) = _ g/35i5s _ (331) 

The rational part 7 3 4,n ° dressing ) cannot be fixed by the principle since the number of inde- 
pendent positively indexed with transcendentality 7 is 64 and we have only 34 spin values 
available. However, with some trial and error, we obtained the following simple formula 

7 (4,no dressing) = I s 7 + 7 S lf6 + 15 S 2 , 5 - 5 S 3 , 4 - 29 5 4)3 - 21 5 5 , 2 - 5 5 fl ,i (3.32) 
-40 5i, liB - 32 5 1<2l 4 + 24 5 1)3 , 3 + 32 5 1|4j2 - 32 5 2)1 , 4 + 20 5 2)2 , 3 
+40 5 2 , 3 , 2 + 4 5 2 , 4 ,i + 24 5 3 ,i, 3 + 44 5 3 , 2 , 2 + 24 5 3)3 ,i + 36 S 4) i, 2 + 36 5 4)2 ,i 
+24 55,1,1 + 80 5i,i,i, 4 - 16 Si,i j3 , 2 + 32 5i,i, 4 ,i - 24 Si, 2)2 , 2 + 16 5i, 2 , 3 ,i 
— 24 5i, 3 ,i, 2 — 24 5i, 3 , 2 ,i — 24 5i, 4 ,i,i — 24 5 2 ,i, 2 , 2 + 16 5 2 ,i, 3 ,i — 24 5 2 , 2 ,i, 2 
— 24 5 2 , 2 , 2 ,i — 24 5 2 , 3 ,i,i — 24 5 3 ,i,i, 2 — 24 5 3 ,i, 2 ,i 
-24 5 3 , 2 ,i,i - 24 5 4 , 1,1,1 - 64 5i,i,i, 3 ,i. 

This formula is minimal in the sense that it is the only solution with coefficients having 
2 as the largest denominator. Notice that there are no solutions with integer coefficients. 
Eqs. < $M 13311 and |3|) are the main result of this paper. 

(n) 

Now, given the above formulas, one can give up the goal of obtaining 73 (s) as exact 
rational numbers. Working with a moderate precision it is possible to go quite further in 

(n) 

the spin. We computed the coefficients 73 ; (s) with 200 digits accuracy up to s = 300. The 
result can be compared with the guessed formulae above. The equality is perfect for all 
digits ! This means that the proposed 73 is actually the only solution compatible with 
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our twist-3 principle. Also, at three loops, one could exclude possible weaker forms of the 
principle admitting nested harmonic sums with some negative index. 

In the next Section, we shall test these expression in the large spin limit showing that 
they reproduce the four-loop universal cusp anomalous dimension. 

3.4 The cusp anomalous dimension 

The cusp anomalous dimension r cusp ((7) is defined in the large s limit of the minimal 
twist- L anomalous dimension as 

7l( s ! 9) = r cus P (5) log s + subleading at s -» 00, (3.33) 
r cU sp(<7) = E r ^ P <? 2n - (3-34) 

n>l 

It is expected to be twist-independent. Its all-loop perturbative expansion is generated by 
the Beisert, Eden, Staudacher equation J4q]. We shall use it as a non-trivial check of our 
finite spin expressions. For completeness, we mention that r cusp (<7) is quite an interesting 
object in the context of AdS/ CFT duality. From the dual string theory, it is known at 



strong coupling at leading and next-to-leading order [£(], |61], [62|]. Checks of this strong 
coupling limit have been recently discussed in [|63|]. 

The large spin expansion of nested harmonic sums can be accomplished by the meth- 
ods summarized in Appendix (A). The calculation is straightforward for the 1, 2, 3 loop 
contributions. Using the results from App. ( |A.4 ), we find 

ifuip = 4, (3.35) 

ifulp = -4S 2 (oo), (3.36) 
Irak = 4 Sf(oo) + 2 5 4 (oo) + 8 5 3 ,i(oo). 

The exact values are 

5 2 (oo) = C 2 = y, 5 4 (oo) = C4 = ^, S3,i(oo) = C4 + ^C 2 2 = ^- (3-37) 
Collecting, we find the correct result 

r cusp (<?) = 4 5 2 - 2 -t- S 4 + ^ 9* + 0(g 8 ). (3.38) 



The expansion of the four loop contribution needs some reshuffling. An alternative, 
equivalent form more suitable to study the large spin expansion turns out to be 

7 (4,no dressing) = £ (5 ^ _ 4 ^ S f + (3.39) 

' 3 

+4(5 5 5 - 10 5 3 , 2 - 145 4 ,i + 16 5 3 ,i,i) S\ + 

+(32 Sf + 8 S 2 Si + 47 S 6 + 88 S 2 5 3 ,i - 64 5 4 , 2 + 48 S 5 ,i - 24 5 2 , 2)2 

-104 5 2 , 3il + 120 5 4 ,i,i " 192S 3ll ,i,i) 5 X 

-2 5 3 5| - S 5 S 2 - 4 5 3)2 S 2 + 60 5 4 ,i S 2 - 120 5 3 ,i,i S 2 

289 189 S 7 256 
— ^- ^3 S 4 —L - — 5 3 5 3 ,i + 136 5 4 , 3 - 24 5 5)2 - 32 5 6 ,i - 64 5 2 , 4 ,i 

+64 5 3)2]2 + 80 5 3)3] i - 80 S^^i + 128 5 2)3) i,i - 128 S^i^i + 256 53,1,1,1,1. 
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In this form, the techniques of Appendix (A) are enough to expand at large s the four loop 
anomalous dimension. Including the dressing term, the result is 

itulp = +4CI - 8PCz- (3.40) 

With the standard choice = £3, we obtain 

r cusp ( 5 ) = 4 5 2 - — g* + — f - + Hi J / + • • • • (3.41) 

It agrees with the analytical prediction by Beisert, Eden and Staudacher J46|] as well as 



with the independent numerical prediction of [ |47| ] from a remarkable four loop Feynman 
diagram computation. 

Notice that Eq. fl3.40D comes from many contributions and requires several cancella- 
tions of higher powers of log s. It tests rather severely our 4 loop proposed expression. 

The next-to-leading contributions are also interesting and shall be discussed in a 
forthcoming publication [|66|]. Here, for completeness, we quote some of them. The ex- 
pansion has the general form 



00 ^ k 

7 = plogs + ^3^ p M log* S, (3.42) 
fc=o e=o 

1 ~" (3.43) 



The first terms are 

? 2 7T 2 4 11 7T 4 fi / 73 7T 6 >2 „, \ fi 

Poo = -2 Cs 3 4 + ^y Cs - Cb j 5 6 + constant • g 8 + • • • , 

2 27T 2 4 / llTT 4 \ 6 / 737T 6 4 7T 2 >2 a> \ o 

"° = 4 » - — » + (-ir - 4& j 9 + (-IN + — 6 +4& - Wa - 2& j 9 + 



Also, we computed 



P22 = -8 5 6 + 47r 2 5 8 + 

32 o 
P33 = y 5 + ■ ■ ■ 



The coefficient of log 5, i.e. p is precisely the cusp anomalous dimension. 
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4. Conclusions 



In summary, we have considered twist-3 operators with spin s in the si(2) sector of M = 4 
SYM built out of three scalar fields with derivatives. We have extracted from the Bethe 
Ansatz equations the exact lowest anomalous dimension 7(5) of scaling fields for several 
values of even s. Here, exact means in analytical form as a known coefficient at each loop, 
up to the four loop level where wrapping problems invalidate the Bethe equations. 

From these results, we have been able to provide closed expressions for the spin de- 
pendence of 7(s) up to the four loop level and including the contributions from the dress- 
ing phase. The expressions satisfy a rather simple and new transcendentality principle ex- 
tending at twist-3 the KLOV idea p7|]. As an application, we have computed the large 
s limit of 7(s) and checked that the four loop universal cusp anomalous dimension is 
reproduced. 

The availability of the full spin dependence allows to test generic features of anoma- 
lous dimensions in conformally invariant planar field theories in the spirit of [ |64| , |65| ]. 
More detailed results on this important topic will appear in a forthcoming publication [Fj6p. 

From the point of view of integrability, it is very interesting to test the consequences 
of the complicated dressing phase Eq. fl2.1ip at finite spin. Indeed, the only weak-coupling 
test of of the conjectured dressing phase has been the calculation of the cusp anomalous 
dimension at infinite s. We consider a relevant result that the leading weak-coupling 
contribution from dressing nicely fits the twist-3 transcendentality principle at finite spin. 
From this point of view, the choice of twist-3 has been crucial to push wrapping effects 
beyond four loops where the dressing first appears. 

As a final comment, we remark that twist-3 operators in QCD built out of quarks and 
gluon fields have a well-known phenomenological relevance in polarized DIS applica- 
tions [0, 58 1 . Here, we have considered the scalar sl(2) sector which is not necessarily 



related to the more QCD-like channels. However, we believe that the ideas presented in 
this paper will also be useful in other twist-3 sectors, possibly exploiting supersymmetry 
to partially connect different channels. 
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A. Harmonic sums 

We collect in this Appendix some useful properties of (nested) Harmonic sums. General 



useful references can be found in [59]. 



-15- 



A.l Definition 

For any integer N and multi-index a = (a\, . . . , a^) £ Z k , we define recursively 

s *w = £ ^. ( A -D 

n=l 

S.xW = E^Sx(n). (A.2) 

n—l 

In the following we shall need only the Oj > case. 
A.2 Some shuffle algebra relations 

The product of a simple S a and a nested sum 5b can be written 

S a Sbi,...,b k = S a ,bi,...,b k + Sb lt a,b 2 ,...,b k + ■ ■ ■ + Sb x ,...,b k ,a (A.3) 
— Sa+b u -,b k — Sb 1 ,a,+b 2 ,...,b k — ■ ■ ■ — Sb lt ... t a+b k - 

In particular 

5a Sb = S a b + 5ft a — S a +b- (A.4) 

These relations can be used to reduce sums of the form S a ... a - One finds for instance 

S aa = \{S 2 a + S 2a ), (A.5) 

Saaa = l{S 3 a + 3 S a S 2a + 2 5 3a ), (A.6) 

Saaaa = ^ + 6 S 2 a S 2a + 3 5 2 2 a + 8 5 a 5 3a + 6 5 4a ) . (A.7) 

A.3 Linear relations 

Given a particular nested sum 5 a , there are linear relations between all the sums 

S a i, a' = 7ra, 7r = permutation, (A.8) 

involving also 5 sums with a smaller number of indices. These linear relations are ob- 
tained from the equations (one for each permutation a') 



S a > S a > 2t a' a ,a' A „„ = shuffle relation Eq. (|A3> (A.9) 
For instance, taking a = (1,1,3) we find 

-Si £l,3 = 51,1,3 + 51,1,3 + £l,3,l — £2,3 — £l,4, (A.10) 

£l £3,1 = £l,3,l + £3,1,1 + £3,1,1 — £4,1 — £3,2, (A.ll) 

£3 £1,1 = £3,1,1 + £1,3,1 + £1,1,3 — £4,1 — £l,4- (A.12) 
They can be used to obtain two of the three 3-index sums in terms of the third. 
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A.4 Evaluation of 5 a at N = oo 

We define the generalized, finite N, (-functions 



W>ni>n 2 >--->n r >0 1 r 

The values of Z at i\T = oo are the multiple zeta functions 

^a(OQ) EE Ca- (A.14) 

Multiple zeta functions can be reduced in terms of (known) elementary zeta functions. 
The relations between Z and S-sums are trivial. For instance 

Sa(oo) = Ca, (A.15) 

S a ,b{°°) = Co,6 + (a+b, (A.16) 
S a ,b,c{°o) = ( a ,b,c + (a+b,c + Ca,6+c + (a+b+c- (A.17) 

The general case is obtained by summing over all possible (a obtained by splitting the 
multi-index of S in order-respecting groups and taking the sum within each group. For 
instance 

5a,&,c,d(o°) = Ca,b,c,d + (a+b,c,d + (a,b+c,d + Ca,b,c+d + Ca+b,c+d + 

+Ca+b+c,d + Ca,6+c+d + (a+b+c+d- (A.18) 

Some examples relevant to compute the three loop cusp anomalous dimension are 
the following. From 

Cs,2 = -yC5 + 3C 2 C3, (A.19) 

£3,1 = ^2, (A20) 

C4,i = Cs.i.i = 2C 5 - C2C3, (A.21) 

we obtain 

53,2(00) = -|Cb + 3C 2 C3, (A.22) 

^3,1(00) = C4 + ^C 2 2 = ^, (A23) 
5 4 ,i(oo) = 3 C5-C2C3, (A.24) 
53,1,1(00) = -|C5 + C2C3- (A.25) 

A.5 Asymptotic expansions of harmonic sums with positive indices 

We first define 



n=l 
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The following expansions hold (Bk are Bernoulli's numbers) 

5 lW =logtf + 7s + ^-£^, (A.27) 

fc>l 

, * a-2iV-l 1 (2A; + a-2)!S 2fc „ T 

W^^-M^W^ 1 a 6 N,a>l.(A.28) 

Taking derivatives with respect to a we obtain immediately expansions for S^\n). 

Multiple (nested) sums S a can be evaluated as follows. Let a = (ai, 02, . . . , a*.). Sup- 
pose that the expansion of 5 02i ... jafc (JV) is known. Its general form will always be 

Sa 2 ,...,a k = ^2 Cp,q J7Z ■ (A.29) 



p,q 



Replacing this expansion in 



N 1 



*a=£-^ ( A - 3 °) 



1 n 

n—l 



we obtain 



log p n 



*a = £<™£^? (A.31) 

P;9 n=l 

Usually this determines the expansion of the sum apart from the constant term. This is 
5 a (oo) and can be evaluated by the methods discussed in the previous sections. Often, 
it is useful to reduce internal sums of the form S aaa ,,, a with the general formulae that we 
also discussed. 

A.6 Counting fixed transcendentality terms 

Let us consider a given index set a = (ai, . . . , a n ) with a t > 0. The various S a / with a' 
being a permutation of a are not all independent due to the shuffle-algebra relations. The 
number of independent sums is counted by the second Witt's formula. Suppose that 

a = (a, ■ ■ ■ ,a,b_-b, ■ ■ ■ ), (A.32) 

and let n = n\ + n 2 H + n^. The number of independent sums with this index set, up 

to permutations, is 

«a) =*.(«,,...,»») = i £ m {ni/d j nW {nN/dr (A33) 

where the Mobius function fi(d) is 

1, d=l, 

fi(d) = 0, d contains the factor p 2 with p prime > 1, (A.34) 
(— l) 5 , d is the product of s distinct primes Pi > 1. 
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Up to transcendentality 7, the cases that occur are (all different letters stand for different 
numbers) 

1(a) = 4(1) = 1, (A.35) 
l( a_-_a ) = l n (ri) = 0. (A.36) 

n>l terms 

(This relation is well known since S aaa ... a can always be expressed in terms of product of 
smaller S) 

l{ a_--_a b) = £ n+1 (l, n) = = 1, (A.37) 

— * — n 1! • n! 

7i terms 



l(aabb) = 4(2,2) = - 



4! 2! 
M1)^T + M(2): 



2! -2! ^ y 1! • 1! 



(A.38) 



l{abc) = 4(1, 1, 1) = ^ TTTTT! = 2 > ( A39 ) 
i(aaa66) = 4(2, 3) = = 2, (A.40) 

^(aaic) = 4(2, 1, 1) = \ 2[ A ^ V = 3- (A.41) 
Now, let us begin with transcendentality 1. There is a single sum with multiplicity I = 1 

a : £ 

1 : 1 

Hence, the number of independent positive sums with transcendentality 1 is 6i = 1. 
At transcendentality 2 we have 

a : I 

11 : 

2 : 1 

Hence, the number of independent positive sums with transcendentality 2 is 62 = 1- 
At transcendentality 3 we have 

a : I 
111 : 

12 : 1 

3 : 1 

Hence, the number of independent positive sums with transcendentality 3 is 63 = 2. From 
the same table we can also count the number iV 3 of (products of) simple sums with total 
transcendentality 3. Let 

(n + k - l\ . A „„, 

Rn,k = k , (A.42) 
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be the number of combinations of n objects in groups of k with repetitions, 
we read 

N 3 = ^1,3 + 6162 + 63 =4. 



From the table, 
(A.43) 



At transcendentality 4 we have 

a : I 

1111 : 

112 : 1 

13 : 1 

4 : 1 

22 : 

Hence, the number of independent positive sums with transcendentality 4 is 64 = 3. 
At transcendentality 5 we have 

a : I 
11111 : 

1112 : 1 

113 : 1 

14 : 1 

5 : 1 
122 : 1 

23 : 1 

Hence, the number of independent positive sums with transcendentality 5 is 65 = 6. From 
the same table we can also count the number JV5 of (products of) simple sums with total 
transcendentality 5. It is 

N 5 = Rb u B + Rb u 3 h + Rb u 2 63 + 61 64 + 65 + 61 R b2<2 + 6 2 63 = 16. (A.44) 
At transcendentality 6 we have 

a : I 
111111 : 
11112 : 1 
1113 : 1 
114 : 1 
15 : 1 
6 : 1 

Hence, the number of independent positive sums with transcendentality 6 is be = 9. 



a 


I 


1122 


1 


123 


2 


24 


1 


222 






(A.45) 
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At transcendentality 7 we have 



a 


I 








a 


1111111 
1111111 


(J 








1123 


111112 










124 


11113 










25 


1114 










133 


115 










34 


16 










223 


7 




1222 


11122 


2 





(A.46) 



Hence, the number of independent positive sums with transcendentality 7 is b-j = 18. The 
same counting as before, taking into account the non trivial Rb 3 ,2 = #2,2 = 3 (i.e. aa, ab 
and bb) gives 

N 7 = 64. (A.47) 
We have computed N k for the even and odd k up to k = 11 and checked that 

N k = 2 k ~K (A.48) 

The same counting can be done including sums with one or more negative indices. The 
number of sums 6^ for k = 1, . . . , 7 is now 

6^ = 2, 3, 8, 18, 48, 116, 312. (A.49) 

Evaluating , the general formula seems to be 

= 2-3 fc ~ 1 . (A.50) 
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